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Abstract. The low-temperature asymptotic expressions for the Casimir interaction between two real metals 
described by Leontovich surface impedance are obtained in the framework of thermal quantum field theory. 
It is shown that the Casimir entropy computed using the impedance of infrared optics vanishes in the limit 
of zero temperature. By contrast, the Casimir entropy computed using the impedance of the Drude model 
attains at zero temperature a positive value which depends on the parameters of a system, i.e., the Nernst 
heat theorem is violated. Thus, the impedance of infrared optics withstands the thermodynamic test, 
whereas the impedance of the Drude model does not. We also perform a phenomenological analysis of the 
thermal Casimir force and of the radiative heat transfer through a vacuum gap between real metal plates. 
The characterization of a metal by means of the Leontovich impedance of the Drude model is shown 
to be inconsistent with experiment at separations of a few hundred nanometers. A modification of the 
impedance of infrared optics is suggested taking into account relaxation processes. The power of radiative 
heat transfer predicted from this impedance is several times less than previous predictions due to difi'erent 
contributions from the transverse electric evanescent waves. The physical meaning of low frequencies in the 
Lifshitz formula is discussed. It is concluded that new measurements of radiative heat transfer are required 
to find out the adequate description of a metal in the theory of electromagnetic fiuctuations. 



1 Introduction 

During the last few years, complicated problems connected 
with the concept of quantum fluctuations generated much 
interest among specialists in gravitation and cosmology, 
dispersion forces, Bose-Einstein condensation, nanotcch- 
nology, radiative heat transfer and related subjects. Van 
der Waals and Casimir forces, which are different kinds 
of dispersion forces, arise from zero-point oscillations of 
the electromagnetic field and thermal photons. They act 
between closely spaced macrobodies, between a micropar- 
ticle and a macrobody or between two microparticles. The 
theory of dispersion forces is based on quantum statistical 
physics. For real bodies at temperature T described by a 
dielectric permittivity depending only on frequency, the 
van der Waals and Casimir forces acting between them 
are calculated in the framework of Lifshitz theory [ll[2l[3] . 
Originally Lifshitz theory was developed using the concept 
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of an oscillating electromagnetic field and the fluctuation- 
dissipation theorem. Later the main equations of this the- 
ory, including the famous Lifshitz formula, were rederived 
in different formalisms [4,5,6,7 and in particular on the 
basis of thermal quantum field theory in the Matsubara 
formulation [S]. The Lifshitz theory was recently used for 
the interpretation of many experiments on the measure- 
ment of the Casimir force [9, 10, 11, 12, 13, 14, 15, 16,17i|18l 
[T9l[20l[2T[ [22,23,24j. in the application of the Casimir and 
van der Waals forces in nanotechnology [25 l l26 l l27 p 28 p 29j . 
in Bose-Einstein condensation 30,31J and also for the de- 
scription of radiative heat transfer between two bodies at 
different temperatures through a vacuum gap [32p33j. 

The application of Lifshitz theory to real metallic bod- 
ies at nonzero temperature has led to controversial results, 
depending on the used model of dielectric permittivity. 
If the boundary bodies are described by the free elec- 
tron plasma model the ensuing thermal correction to the 
Casimir force [34,35 is in qualitative agreement with that 
obtained for ideal metals on the basis of thermal quantum 
field theory in Matsubara formulation [81136]. If, however, 
metal boundaries are described by the Drude model which 
takes relaxation into account, the thermal correction at 
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short separations is many hundred times larger than for 
ideal metals and two times smaller than the latter at large 
separations [37, 38^. In the case of perfect crystal lattices 
with no impurities, the entropy of a fluctuating field (i.e., 
the Casimir entropy) calculated using the Drude model 
takes a negative value when the temperature vanishes, 
i.e., the Nernst heat theorem is violated ^39,40]. In [HI 
I42j it was argued that for real metals with impurities this 
violation does not occur, but in the case of perfect crys- 
tal lattices the problem remains unsolved. If the dielectric 
permittivity of the plasma model is chosen the Nernst heat 
theorem is satisfied 39,40 . Importantly, the application 
of the Drude dielectric function was found to be inconsis- 
tent with experiment pJl[T7lfT8lfT9ll20] . On the contrary, 
the plasma model approach is consistent with experimen- 
tal data [15] . 

Interestingly, similar problem arises in the theory of 
the thermal Casimir force acting between dielectrics. If the 
static dielectric permittivity of dielectric materials is sup- 
posed to be finite, Lifshitz theory is found [44[l45] to be in 
agreement with thermodynamics. If, however, the dc con- 
ductivity of dielectric materials is taken into account, the 
Nernst heat theorem for the entropy of a fluctuating field 
is violated [iilHS] . The same is true in the metal-dielectric 
configuration [46] depending on the finiteness of the static 
permittivity of a dielectric plate [47ll48]. Thus, thermo- 
dynamics provides a test for the validity of various mod- 
els of material properties: only the thermodynamically- 
consistent models should be used. In this connection it is 
notable that just the Drude model of a metal, which was 
shown to imply a violation of the Nernst heat theorem 
for the Casimir entropy, was found to be in contradiction 
with experiments (controversial opinions on this subject 
can be found in [4^jl50]). 

As an alternative to the dielectric model, the optical 
properties of a metal surface can be characterized in terms 
of the Leontovich surface impedance, together with the 
corresponding boundary conditions 51.. In the framework 
of the Lifshitz theory it was employed in [52j at zero tem- 
perature and in [55[f51] at nonzero temperature. It should 
be kept in mind that both models of real metals, the one 
based on the frequency-dependent dielectric permittivity 
as well as the one using Leontovich impedance, are ap- 
proximations, each having its own range of validity. The 
concept of a frequency-dependent permittivity is inappli- 
cable in the frequency region of the anomalous skin effect, 
where spatial dispersion contributes critically. Physically 
this is explained by the fact that for these frequencies 
the penetration depth of the electromagnetic field inside 
a metal becomes of the same order as the mean free path 
of conduction electrons and remains much less than the 
distance traveled by an electron during the period of the 
field. As a result, the spatial non-uniformity of the field 
renders impossible a macroscopic description in terms of 
a dielectric permittivity depending only on the frequency 
[51j . On the other hand, impedance boundary conditions 
retain their validity at the frequencies of the anomalous 
skin effect because the field inside a metal near the surface 



can be considered as a plane wave propagating perpendic- 
ular to the surface. 

However, Leontovich impedance, though applicable to 
the description of the anomalous skin effect, cannot be 
used at short separations where its magnitude is not much 
less than unity. The physical reason for this is that with 
decreasing separations the relevant characteristic frequen- 
cies enter the optical region, where the field inside a metal 
cannot be considered anymore as propagating perpendic- 
ularly to the metal surface. In this frequency region the 
concept of dielectric permittivity depending only on the 
frequency is satisfactory (because the distance traveled 
by an electron during the period of the field is much less 
than the penetration depth). However, the magnitude of 
the dielectric permittivity is not large enough, and, as a 
consequence, the angle of refraction depends on the an- 
gle of incidence. Thus, the impedance boundary condition 
does not apply. 

The explicit analytic forms of the impedance function 
are available in the asymptotic regions of the normal and 
anomalous skin effect and in infrared optics. As was no- 
ticed without a detailed proof in [54 , the Casimir entropy 
calculated using the impedance of infrared optics van- 
ishes when the temperature goes to zero, i.e., the Nernst 
heat theorem is satisfied. The same is proved in |55) for 
the Casimir entropy calculated with the impedance of the 
anomalous skin effect. We stress that at large separations 
the impedance approach leads to magnitudes of the ther- 
mal Casimir force in qualitative agreement with the case 
of ideal metals. 

A critical problem of the impedance approach is the 
choice of the functional dependence of the impedance func- 
tion on the frequency. In ^54j it was argued that the impe- 
dance function valid in the region around the characteris- 
tic frequency should be extrapolated to lower frequencies, 
all the way to the zero Matsubara frequency. In I56j quan- 
titative arguments were adduced in favor of the statement 
that the use of different impedance functions within differ- 
ent frequency regions in accordance with their applicabil- 
ity conditions would be thermodynamically inconsistent. 
In the first part of this paper we apply the thermodynamic 
test to the impedance function of [33] which provides a 
smooth analytic interpolation between the impedances of 
the normal skin effect and infrared optics. The substitu- 
tion of the interpolated impedance [33] in the Lifshitz for- 
mula presents an explicit example of the situation where 
different impedances are used at Matsubara frequencies 
belonging to different frequency regions. We present a rig- 
orous analytic proof that the entropy of a fluctuating field 
in this situation goes to a positive value when the tem- 
perature vanishes. In other words, the Lifshitz formula 
combined with the interpolated impedance is thermody- 
namically inconsistent. 

The problem of the thermal correction to the Casimir 
force was further discussed in [571 1581159] where the impedance 
of the normal skin effect was used in the computations. In 
this approach at a separation of 1 fim the thermal cor- 
rection was found to be about 30 times larger than for 
ideal metal plates. This was explained by the dominant 
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contribution of the transverse electric evanescent waves 
of rather low frequencies. The same impedance function 
was applied [33^ to compute the radiative heat transfer 
through a vacuum gap between two metal surfaces at dif- 
ferent temperatures. This problem was considered previ- 
ously in [601I6T] (and recently in [32]) using the formalism 
of dielectric permittivity. As proposed in [33], the prob- 
lems of the thermal Casimir force and the radiative heat 
transfer through a vacuum gap are closely related and 
their simultaneous investigation can be very fruitful and 
elucidating. 

In the second part of this paper we consider both the 
thermal correction to the Casimir force and the power of 
radiative heat transfer using different impedance functions 
and compare the results with those obtained previously 
using the formalism of dielectric permittivity and with 
experiment. We try to address, from a purely phenomeno- 
logical point of view, the radiative heat transfer through 
a vacuum gap in connection with the problem of thermal 
Casimir force between real metals. The aim of our work on 
this subject is to come up with a phenomenological model 
for a real metal at room temperature, that takes dissi- 
pation into account, and is at the same time consistent 
with available experimental facts. In view of the existing 
controversies among theoreticians about the correct way 
to do this, it is our view that finding a reasonable empiri- 
cal model would be a valuable guide for further theoretical 
studies. We perform computations in the framework of two 
different formulations of Lifshitz theory along the real and 
the imaginary frequency axis. This permits to specify the 
comparative role of the travelling and evanescent waves in 
the physical phenomena under consideration and to deter- 
mine the application range of different approximations. It 
is shown that use of the impedance of the normal skin ef- 
fect results in enormously large thermal corrections to the 
Casimir force at separations of about 200-300 nm which 
are inconsistent with already performed experiments. As 
an alternative, a phenomenological generalized impedance 
of infrared optics is constructed taking into account relax- 
ation processes specific for this frequency region which are 
not connected with the electron-phonon interactions. The 
thermal correction to the Casimir force computed with 
the generalized impedance function is shown to be quali- 
tatively the same as is known for ideal metals from thermal 
quantum field theory in Matsubara formulation. It is con- 
sistent with all available experimental data. The suggested 
impedance function is applied to calculate the power of ra- 
diative heat transfer between Au plates at different tem- 
peratures. At short separations between the plates the ob- 
tained power of heat transfer per unit area is several times 
less than the one predicted previously in literature using 
the dielectric function or the impedance characteristic for 
the region of the normal skin effect. Thus we find that, 
depending on the chosen model of the metal, one obtains 
largely different predictions for both the thermal Casimir 
force and the power of radiative heat transfer. This result 
underlines the crucial role of new experiments on the pre- 
cision measurements of the Casimir force and the power 
of radiative heat transfer. 



The paper is organized as follows. In Sec. 2 we present 
two equivalent forms of Lifshitz formula, which involve 
evaluating the reflection coefficients for imaginary and real 
frequencies, respectively. In Sec. 3 we give the explicit ex- 
pressions of the reffection coefficients in terms of a fre- 
quency dependent dielectric permittivity and in terms of 
Leontovich surface impedance. In Sec. 4 we demonstrate 
that Lifshitz formula combined with the impedance func- 
tion of the infrared optics is thcrmodynamically consis- 
tent. In Sec. 5 we prove that the substitution of the inter- 
polated impedance into the Lifshitz formula results in a 
violation of the Nernst heat theorem. Sec. 6 is devoted to 
the computation of the thermal correction to the Casimir 
force using different impedance functions. The results are 
compared with the case of plates made of ideal metal and 
with experiment. In Sec. 7 the general expression for the 
power per unit area of radiative heat transfer in terms of 
the surface impedance is presented. Sec. 8 contains the 
computation results for radiative heat transfer with dif- 
ferent impedances and dielectric permittivities including 
the new prediction to be tested experimentally. In Sec. 9 
the reader will find our conclusions and discussion. 

2 Lifshitz formula along the imaginary and 
real frequency axis 

In this section we consider two thick dissimilar plane par- 
allel plates (semispaces) in thermal equilibrium at equal 
temperature T, separated by an empty gap of width a. 
Let the z axis be perpendicular to the plates. The Lif- 
shitz formula [T] represents the van der Waals and Casimir 
free energy and force per unit area (i.e., the pressure), 
acting between the plates marked by the indices (f ) and 

f 1 2) 

(2), in terms of the refiection coefficients r^-^' {to, k±) and 

r^g {cj,k±) for two independent polarizations of electro- 
magnetic field {lo is the frequency and k± is the mag- 
nitude of the projection of the wave vector in the plane 
of the plates) . The transverse magnetic polarization (TM) 
means that the magnetic field is perpendicular to the plane 
formed by k± and the z axis, while for the transverse 
electric polarization (TE) the electric field is perpendicu- 
lar to this plane. As was mentioned in the Introduction, 
in the literature there are many different derivations of 
the Lifshitz formula in the framework of quantum statis- 
tical physics, thermal quantum field theory in Matsubara 
formulation and scattering theory (see, e.g., [4l|5l[6l[7l[8] ) . 
The final results of all derivations are represented in one of 
two different forms, as a summation over the Matsubara 
frequencies along the imaginary frequency axis or, alterna- 
tively, as an integral over real frequencies. We begin with 
the more often used representation in terms of imaginary 
frequencies, where the Casimir free energy and pressure 
are given by 

Ha, T) = ^ ^ f - -6i, / k^dk^ (1) 

a=TE,TM 
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P(a,T) = - 



1=0 



10 



E 

c«=TE,TM 



2 

2aqi 



kj_dkj_qi 



(2) 



Here ^; = 2TrkBTl/ti are the Matsubara frequencies, fcs is 
Boltzmann constant, I is a non- negative integer and 



qiuj,kj_) = \jk\- — , 



(3) 



qi EE qii^i,k±_) = \ k\ + 



(1 2") 

The exphcit form of the reflection coefficients r-Q is dis- 
cussed in the next section. Representation ([2]) is conve- 
nient in numerical computations due to fast convergence 
of both the sum and the integrals. 

Using the Abcl-Plana formula [5115^ 

f2{^-\5K^m = F{t)dt (4) 

F{it) - F{~it) 



1=0 



+ i I dt- 

10 



^2^t _ I 

where F{z) is an analytic function in the right half-plane, 
Eqs. ([T]) and ([2]) can be identically rearranged in the form 

T{a,T) = E{a) + AT{a,T), 

P{a,T)^Po{a) + AP{a,T). (5) 

Here £'(0) is given by 

E{a) = ^ / dn k^dk^ (6) 



Q=TE,TM 

with q = k\^) defined in Eqs. The other contribu- 
tion on the right-hand side of the first equality in Eq. ^ 
can be represented in the form 

iksT 



AT{a,T) = 



27r 



(7) 

e — 1 



where <?(x) = (pTuix) + (1>te{x) and 



TM.TE 



(x) 



k \ dk \ 



(8) 



X In 



1 ''tM Te(*'''' ^J-)''' 



(2) (■ h \ 
TM.TEV*'''' Ki^je 



2a yi^^ 



In a similar way, the quantities Po{a) and AP{a, T) in the 
second equality in Eq. ^ are equal to 

H -1 



Po{a) 



' d^ / k±dk±q 
Jo 



E 



„2aq 



Q=TE,TM L 

AP{a,T). 



ri^\i^,k±)r^^\i^,k^) 



1 



(9) 



iksT 



p27rt 



where F{x) = i^TM(a;) -f i^TE(a;) and 

-FTM,TE(a;) = / kj_dk±q\ k\ + — 



(11) 



(2) 



- f 



Note that the quantities E{a) and Po{a) in Eqs. ([6]) and 
© are often called in literature the Casimir energy and 
pressure at zero temperature, and Z\^(a, T) and AP(a, T) 
in Eqs. ([Tj) and (fTO|) are referred to as the thermal correc- 
tions to them. This terminology is, however, correct only 
for plate materials with temperature independent prop- 
erties. In this case the Casimir free energy and pressure 
depend on temperature only through the Matsubara fre- 
quencies and the thermal corrections defined as T{a, T) — 
T{a, 0) and P(a, T) - P(a, 0) coincide with AT{a, T) and 
Z\P(a,r) in Eqs. ^ and If, however, the proper- 

ties of a medium (for instance, the dielectric permittiv- 
ity) depend on the temperature, then the thermal cor- 
rections J^{a,T) - jr(a,0) and P(a,T) - P(a,0) do not 
coincide with A!F{a,T) and AP{a,T). Even in this case 
Eq. (O can be used to compute the total Casimir free en- 
ergy and pressure. In so doing, the quantities E{a) and 
Po(a) in Eqs. (O and ^ can be interpreted as the con- 
tributions to the total Casimir free energy and pressure 
due to zero-point oscillations [they may depend on the 
temperature as a parameter, i.e., in fact E{a) = E{a,T) 
and Po(a) — Po{a,T)] and the quantities AT{a,T) and 
AP{a, T) in Eqs. ([7]) and (fTU]) as the contributions from 
thermal photons. 

It is useful to express the Casimir free energy and pres- 
sure as the integrals over real frequencies lo: 



Ha,T)^-l d.,fc,coth(^) 



(12) 



a^TE.TM 



oc poo 



P(a,T).-^y^ d.j^ d.,fc,coth(J^)(13) 



X Re < fcz 

a=TE,TM 



,-2ife. 



1 - 



where 



K_iiJ, k^) = ^LO^/c-^ -kl= iq{u, k^) . (14) 



We note that real values of k^ correspond to propagat- 
ing waves (PW), while imaginary values of k^ describe 
evanescent waves (EW). Upon using the identity 



coth(a;/2) = 1 -|- 



exp(a;) — 1 



(15) 



we see that the formula for the Casimir pressure can be 
expressed as the sum of two terms, like in Eq. ([5|) 



P(a,T) = Po(a,T) + ZiP(a,r) 



(16) 
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where 



Po{a,T)^-—j diuj dk^k^ 



(17) 



Rc\k, 

Q=TE,TM 



H -1 



1 - 



ra\uj, k±)ra'{uj, kj_) 



(2), 



and 



h f°° f°° 1 

AP{a.T)^ dw dkj_k±_ r 

-'Jo Jo exp(^)-l 



a=TE,TM . 



1 - 



ra\uj, kj_)ra'iuj, k±) 



(2), 



(18) 



Similar representations can be easily obtained for E{a, T) 
and A!F{a, T). As was already noted above, PQ{a, T) phys- 
ically represents the contribution to the pressure from 
zero-point fluctuations of the electromagnetic field, and in 
general it depends on the temperature, because the per- 
mittivities or surface impedances of the plates are tempe- 
rature-dependent. As for AP{a,T), it represents the con- 
tribution from thermally excited electromagnetic fields, 
and it vanishes for T = 0. For this reason, in what fol- 
lows we shall conventionally refer to AP{a,T) as to the 
thermal correction to the Casimir pressure. We further 
consider the decomposition 



AP{a, T) = APpwia, T) + AP^^{a, T) 



(19) 



where Z\Ppw(a,?') and Z\PEw(a, T) represent the contri- 
butions from PW and EW, respectively. It is easily seen 
from Eq. that Z\Ppw(a,T') can be written as: 

^-Ppw(a,r) = / dw r / dk^kl 

" J' exp (^) - 1 Jo 



E 

Q=TE,TM 



Re 



exp 

2ikz o. 



(2), 



(20) 



For EW it holds k± > uj/c and Im(fc2) = g. As a result 
APew{o-, T) takes the form 



APEw{a,T)^— / 
Jo 

X ^ Im 



duj- 



1 



exp 



=TE,TM 



\kBT J 



„2qa 



- 1 Jo 



dqq 



r^\u;,k±)r^' {uj,k±^) 



(2), 



(21) 



It should be noted that APEwiO'jT) vanishes in the case 
of ideal metals, because for 

^ ^ — 1 tbe quantity be- 
tween square brackets in the above equation is real. In 
fact, more generally APew{<i,T) vanishes whenever the 

product of the reflection coefficients is real and less 
or equal to unity, because then the quantity between the 
square brackets in Eq. (j21|) is real and has no zeroes. This 
is the case for example for TE EW, if the metal plates are 
described by the plasma model (see below). 



3 Reflection coefficients in terms of the 
dielectric permittivity and the surface 
impedance 

In Lifshitz theory the material media are described by di- 
electric permittivities that depend only on frequency [Hill 
13] • The description of the dielectric properties of a medium 
by s{uj) takes full account of temporal dispersion but ne- 
glects possible contributions to the van der Waals and 
Casimir force from spatial dispersion. In the formalism of 
the imaginary frequency axis [see Eqs. ((T]) and (O] the re- 
flection coefficients are expressed in terms of the dielectric 
permittivity as follows: 



(") 



(22) 



where n = 1, 2 for the first and the second plates, respec- 
tively, and 



kl 



e(«)(tj)- 



(23) 



fc(")^fc(")(ze,,fci)-^A:i+£(")(zeo|. 

In the formalism of real frequency axis used in Eqs. (fT2 
and the reflection coefficients are just (|22p but calcu- 
lated at real frequencies 



(") 
''tm 



?'te(^'^-l) = 



{uj, k±_) 



£<")(^)g(u^,fc^)-fc^"n^.fc-L) 

£(")(cj)q(tj, k^) + k("\uj, kx_) ' 
fc(")(a;,fc^)-g(c^,fc^) 
fc(«)(w,fci) -f 9(w,fci) ■ 



(24) 



The reflection properties of electromagnetic waves on 
metal surfaces are often described in terms of the Leon- 
tovich surface impedance '51J . For isotropic metal surfaces 
the surface impedance relates the tangential components 
of electric field and magnetic induction in the same way 
as in a plane wave propagating in the interior of a metal 
perpendicular to its surface [5l 



Et = Z{uj) [Bt X n] 



(25) 



Here n is the unit vector normal to the surface and di- 
rected inside the medium. For an ideal metal Z = Q. The 
boundary condition (|25p is valid when \Z\ ^ 1. For good 
conductors this inequality is satisfied within a wide fre- 
quency region. Equation (|25p permits to determine the 
electromagnetic field outside the metal, without consider- 
ing the propagation of electromagnetic waves in the metal 
interior. It is close in spirit to the original Casimir ap- 
proach for ideal metals and to so-called "nonlocal" bound- 
ary condition |63| implied by a dielectric permittivity de- 
pending only on the frequency. The latter condition is in 
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fact equivalent to the standard continuity boundary con- 
ditions in classical electrodynamics but does not require 
to consider field propagation inside a metal. Both the 
standard continuity conditions and the "nonlocal" bound- 
ary condition are based on the spatially local relation 
D{x,oj) = e{uj)E{x,uj) which assumes space homogene- 
ity. Therefore these boundary conditions do not take into 
account the effects of spatial dispersion. As a consequence, 
the standard Lifshitz formula is applicable only in the ab- 
sence of spatial dispersion. An advantage of the surface 
impedance, as compared with dielectric permittivity, is 
that it permits |52j to apply the Lifshitz formula in the 
region of the anomalous skin effect, where the spatial ho- 
mogeneity is violated and the effects of spatial dispersion 
should be taken into consideration. In '6^ the applicabil- 
ity of the condition (P5)) in the region of the anomalous 
skin effect is demonstrated from the solution of kinetic 
equations. In this frequency region, a metal cannot be 
characterized by a dielectric permittivity depending only 
on frequency and, thus, the standard Lifshitz formula is 
not applicable [65] (the Leontovich impedance, however, 
cannot be used at short separations between the plates 
where the condition \Z\ ^ 1 is violated, see Introduc- 
tion). In the frequency regions where both quantities e{uj) 
and Z{uj) are well defined it holds 



are 



Z(^) = 



(26) 



In optics of metals the reflection coefficients are usu- 
ally expressed in terms of Z(uj) rather than s{(jj) [51j . In 
[5^ the Leontovich surface impedance was used to express 
the reflection coefficients in the Lifshitz formula at zero 
temperature. The thermal Casimir force was presented in 
terms of the surface impedance in [S3j. The derivation of 
the Lifshitz formula starting from impedance boundary 
condition ([25| is contained in [54]. Importantly, as was 
shown in 54J , the use of the Leontovich impedance leads to 
different results for the thermal Casimir force between real 
metals, than are obtained by the use of the Drude dielec- 
tric function in [37' . The thermal Casimir force computed 
within the impedance approach was demonstrated to be 
in qualitative agreement with the case of ideal metals and 
in accordance with the fundamentals of thermodynamics 
and with experiment. (The controversies between different 
approaches to the thermal Casimir force are discussed in 
detail in [HSOlIMlEZl ) ■ 

Using the Leontovich surface impedance instead of the 
dielectric permittivity, the reflection coefRcients in Eqs. ([T]) 
and ([2]) (in the formalism of the imaginary frequency axis) 
are given by 



?'te(^^''^-l) 



Cg;-Z(")(ze06 

cg, + Z(")(iC06' 
e,+c(7,Z(»)(zeO' 



(27) 



In the formalism of real frequency axis, the reflection 
coefficients expressed in terms of the Leontovich impedance 



J") 
'tm 



(«) 

'^TE 



, , cfc4c^,fc^)-Z(")(cj)co 

^'^'^^ cfc,(a>,fci) + Z(")(t^)a.' 

, , UJ-ck,{LU,kl_)Z^''\uj) 

^'^'^^ c^ + cfc,(c^,fci)Z(")(c^)' 



(28) 



where kz is defined in Eq. (|14p . An important problem 
arising in the case of real metals is the adequate choice of 
the functions e(w) and Z{ijj). 

The calculation of the surface impedance over the whole 
frequency axis is based on kinetic theory [64! . Here we are 
interested in two frequency regions. One of them is the 
infrared optics defined by the inequalities 



^«(5, «;(r), 

UJ 



UJ < Wp, 



(29) 



where is the Fermi velocity, 1{T) is the mean free path 
of a conduction electron, 5i = c/ujp is the skin depth and 
ujp = 2'Kc/\p is the plasma frequency. The other frequency 
region of our interest is the region of the normal skin effect 
characterized by the inequalities 



;(r) <^Ar(c^,T), 



m « ^, 

UJ 



(30) 



where Sn{uj,T) = c/^27r(To(T)a; and cro(r) is the static 
electric conductivity. 

In the frequency regions ([29]) and ([30]) , both the dielec- 
tric permittivity and the impedance have definite physical 
meanings and are connected by Eq. ([26]). In the region of 
the infrared optics we have 



ei{uj) = 1 - 



UJ 



(31) 



For the normal skin effect we have 



ejv(w) = eAr(cj,T) = i 



47rao(r) 



Zn{uj) = Zn{uj, T) = (1 - i) 



87rt7o(T) 



(32) 



Let us consider two plates at room temperature T = 
300 K at the same separation distance of a few hundred 
nanometers, as in [33]. In this case, the respective char- 
acteristic frequency of the Casimir effect f2c — c/(2a) be- 
longs to the region of infrared optics ([29]l . (As an exam- 
ple, at a = 400 nm, He — 3.75 x lO^'^rad/s and for gold 
ujp « 1.37 X lO^^rad/s and vp « 1.78 x lO^m/s.) In the 
Lifshitz formulas ([1]), ^ and the frequencies of 

order ujp and higher practically do not contribute to the 
result. AH Matsubara frequencies ^/ in Eqs. ^ and ([2]), 
which contribute to the result essentially, belong to the 
region of infrared optics (|29)) with the exception of = 0, 
which belongs to the region of the normal skin effect (|30]) 
(recall that = 2.47 x 10^''rad/s). Real frequencies in 
Eqs. (fT2| and ([13]) . contributing to the result, also belong 
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to the region of infrared optics ((29|) and to the normal skin 
effect CT - 

In j33j the correlation functions of a fluctuating elec- 
tromagnetic field were expressed in terms of the surface 
impedance. The developed formalism was applied to de- 
rive the Lifshitz formulas (fT2|) and (|13|) along the real 
frequency axis and to describe the radiation heat trans- 
fer between two semispaces at different temperatures. We 
consider the impedance function defined by the Drude di- 
electric function which was used in 1331 



1 



U![ll! + i^(T)] 



(33) 



Here 7(T) is the relaxation parameter connected with the 
above used parameters by the equation [68] 



4wo(r) 
r{T) ■ 



(34) 



where r = I/7 is the relaxation time. In the region of 
infrared optics we have 7(1^) <^ As a result, Zu and 
£u in Eq. fSS]) coincide with Zi and £i in Eq. ([3T|) . re- 
spectively. At small frequencies, on the contrary, one can 
neglect w, as compared to 7, and from Eq. (|34|) Z^ and eo 
in Eq. (l33|) coincide with Z„ and £„ in Eq. ((32|) describing 
the frequency region of the normal skin effect. Therefore, 
Eq. (j33p provides an expression of the impedance which is 
valid in both frequency regions of the normal skin effect 
and infrared optics, and it represents a smooth analytic 
interpolation between the two regions. 

Below we demonstrate that the impedance of infrared 
optics, Zi, extrapolated to all lower frequencies, includ- 
ing zero frequency, leads to zero entropy of a fluctuating 
field at T = (Sec. 4). At the same time, the entropy 
of a fluctuating field calculated using the impedance of 
the Drude model in Eq. (|33p approaches a nonzero posi- 
tive value when the temperature vanishes (Sec. 5), hence 
violating thermodynamics. 



4 Thermodynamic test for the surface 
impedance of infrared optics 

We consider the free energy of a fluctuating field given by 
Eqs. H]) and with reflection coefficients ^T^j, (US]) and 
the impedance function (j3ip . Our aim is to find the asymp- 
totic behavior of the free energy and entropy of a fluctuat- 
ing field at low temperatures at separation distances be- 
tween two similar plates of a few hundred nanometers, so 
that the characteristic frequency J7c belongs to the region 
of infrared optics. The perturbation expansions in powers 
of the small parameter k = AnksaT/ (he) can be conve- 
niently carried out by using the dimensionless variables 



kI, y = 2aqi 



(35) 



In terms of these variables the free energy of the fluctuat- 
ing fleld, Eq. takes the form 

1=0 ^ ^ 

X J2 ln[l-r^(zO,2;)e-^]. 

a=TM,TE 

Using the variables (I35p . the reflection coefficients I^T^j are 

y - Z{iQnc)Ci 



rTMiiChV) 



y + ZiiQficKi 

Ci - yZ{iCif2c) 



(37) 



Ci + yZ{iCin,) ■ 

Here, the impedance function of infrared optics (pij) is 
given by 

pCi 



(38) 



where p = Ap/(47ra) = di/{2a) is much less than unity 
throughout the entire region of application of the impedance 
approach. 

In terms of new variables the Casimir energy at T = 
in Eq. ([6]) is given by 



Eia) = 



he 



327r2a3 



00 />00 



JC 



dC / fiC.y)dy (39) 



and the thermal correction to it ([7]) by 

, ^, ihcK f°° , <l>(iKt) — <l^(—iKt) 



327r2a3 ^ 62^^* - 1 

where the following notations are used: 

/(C,2;) = yln [l-4^{t(:,y)e-y 
+ y\n [l-4E«,y)e"1 , 



(41) 



^(x) = / dyf{x,y). 



Now, we expand the function / defined in Eq. (|4T|) in 
powers of the small parameter p using Eqs. ([37| and ((38|) : 



f{x,y) = 2y\n{l^e-y)+4p 



X^ + ?/2 

ev -1 



y[ey - 1)2 



(42) 



Substituting Eq. (gH) into the definition of <1> in Eq. (gl 
we obtain 



$ix) = /o(a:) -I- phix) + p^hix) + 0(p3), (43) 



where 



poo 

Io{x) = 2 ydy\n{l-e-y), 

J X 

r°° x2 + ?;2 

h{x)=i dy^ 



ev -1 



(44) 
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The integrals Iq (x) and Ii (x) are easily calculated 
and are given by: 

/o(a;) = -2 [Li3 (e^) + xLia (e^)] 



= -2C(3) + ^x^ ~ 



1 



\nx + ^x^ + 0{x'^), (45) 

I - 

]-x^ \xix + \x^ + 0{x^) 



h{x) = 8 [Lis (e-^) + xLi2 (e^^) - x^ In (l - e"")] 



C(3)-T^ 



where Li„(z) is the polylogarithm function and C(a;) is 
the Riemann zeta function. As to the integral 12(2^) in 
Eq. pi)l . in is easily seen that 



/2(x) = -48C(3) + 0(a;^) 



(46) 



and, thus, I^ix) does not contribute to 'P{iKt) — <P{—iKt) 
in the perturbation orders under consideration. 
Substituting Eq. (|45]) into Eq. ((431), we obtain 



The asymptotic behavior of the entropy of a fluctuat- 
ing field defined as 



dT{a,T) 
df 



S{a,T) 



can be found by differentiating Eq. ([49 



(50) 



S^{a,T) 



3kB 



- 2C(3) 
a 



T 

Toff 
167r3 



\ . 47r3 T 

135 ieff 



T 



135 Teff 



(51) 



As is seen in Eq. (I5ip . entropy goes to zero when temper- 
ature vanishes in accordance to the Ncrnst heat theorem. 
Thus, the Leontovich impedance of the infrared optics 
withstands the thermodynamic test. The Lifshitz formu- 
las ([H and p2)) combined with the impedance of infrared 
optics are shown to be consistent with the requirements 
of thermodynamics. 



<l>{iKt) - <l>{-iKt) = ni^Kt)'^ ~ - i{Kt)^ 



(47) 



+ 4:P 



3 



O [{ntY 



Now, it is easy to calculate the free energy of the fluc- 
tuating field at smaU k from Eqs. ©, (gO]) and (gT]) 



T,{a,T)^E{a) - 



720a3 



45C(3) 
87r6 ' 



1 



167r4 



45p fm 3 

'^4 1^27r2'' 



90 



(48) 



It is convenient to introduce the so-called effective tem- 
perature ksTcs = tific and to use the effective penetra- 
tion depth of electromagnetic oscillations into a metal (|29p 
in the frequency region of infrared optics, Si = c/ujp = 
Ap/(27r). Then, Eq. (|48p can be rearranged in the form 



T{a,T) = E{a) 



TT^hc /45C(3) / T 



720a3 I 7r3 \T,s J \T,s 



T 



a 



90C(3) f T 



T 

Trff 



(49) 



In the above, we have restricted our consideration to the 
second perturbation order in the small parameter p. In the 
same way as in [70] it can be shown that the higher pertur- 
bation orders in p contain only terms of order 0(k") with 
n > 5. It is notable that Eq. ([35]) coincides [MlliO] with 
the free energy of the fluctuating field at low temperatures 
obtained from the Lifshitz formula combined with the di- 
electric permittivity of the plasma model £^(0;) in Eq. ([?T]) . 
Thus, the characterization of a metal by means of the di- 
electric permittivity and Leontovich surface impedance in 
the frequency region of infrared optics leads to the same 
asymptotic behavior of the free energy at low tempera- 
tures. 



5 Thermodynamic test for the surface 
impedance of the Drude model 

In the previous section we have extrapolated the impedan- 
ce function of infrared optics to all lower frequencies, in- 
cluding zero frequency. However, in the Lifshitz formula, 
Eq. ([1]), the zero Matsubara frequency is situated out- 
side the region of infrared optics. In a similar manner, 
small real frequencies in the Lifshitz formula ([T^ sat- 
isfy inequalities ([50)1 and, thus, belong to the region of 
the normal skin effect. Because of this, it seems reason- 
able to use the impedance function (|33p which coincides 
with the impedances of infrared optics and of normal skin 
effect at high and low frequencies, respectively, and pro- 
vides a smooth analytic interpolation between the two fre- 
quency regions. In this section we put the Lifshitz formu- 
las, Eq. ^ and Eq. ([T^ . combined with the impedance 
of the Drude model ([33]) to a thermodynamic test. 

For I > 1 one can introduce the dimensionless variable 
xi — ^(T)/^i and rearrange the impedance of the Drude 
model ([55]) in the form 



pCi 



1 

l+xi 



(52) 



If the relaxation parameter 7(T) is equal to zero, = 
and Eq. ([5^ coincides with Eq. ([55]) for the impedance of 
the infrared optics. It is easily seen that for metals with 
perfect crystal lattices a;; ^ 1 at sufficiently low T. In 
fact, for T = 300K for good metals it holds 7 - 10" - 
lO^'' rad/s (as an example, for gold 7 = 5.32 x 10^^ rad/s), 
whereas = ^il and ^1 = 2iTkBT/h = 2.46 x lO" rad/s 
leading to xi — 7/^; < 0.22. When T decreases from 
T = 300 K to approximately Td/A, where To is the Debye 
temperature (for gold Td = 165 K [Tlj), 7(T) decreases 
linearly, ^{T) ~ T, i.e., following the same law as ^j. At 
T < Td/A the relaxation parameter decreases according 
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to the Bloch-Griineisen law, ^{T) ~ T^, due to electron- (f38| . with omitted contribution from the zero Matsubara 
phonon collisions ^68j and as j{T) ^ at liquid helium frequency. As a next step, we expand the integrand in 



temperatures due to electron-electron scattering [71] ■ At Eq. ((58|) in powers of a small impedance of infrared op- 

T = 30K and lOK it holds j{T)/^i{T) x 4.9 x lO'^ and tics, ZiiiQn^), and retain only the zero order contribution 

1.8 X 10~^, respectively. The magnitude of parameter xi [recall that Zi{iC,ifIc) goes to zero when T vanishes]: 
decreases further to zero with T 0. 

We represent the free energy of the fluctuating field in J^^^"^^ (a, T) ~ T-''^^^ (a, T) 

Eq. in the form , „ oo ^oo , / 2 

fesT^, r dy ( (■ , y 

Tn{a,T)=Tt'\a,T)+T^^'\a,n (53) 2™^ ^ 7,, 1 - 



where we separate the terms with zero and nonzero Mat- 



0[Z,{iCif2,)ii,if]. (59) 



subamfrequencies. Substituting the impedai^ function, 3^^^. definition of £1 in Eq. m, the sum 

Eq. m the refiection coefficients, Eq. one ob- ^ ^^^^^ ^^^^ 



tams 

„2 /n „.\ „2 



^TM(0,2/)=r-TE(0,y) = l, (54) 7(r) 2 2^-3/2 T dy ( y^ 



which leads to ^ - ^ 2^^^ + P > J + Z 



1=1 



^{1=0) f _ ksT , _y. (60) 

•^D i - g^^2 y^-y^'^K^ e ) . [00) expand the integrand in Eq. ^ in powers of 

and introduce the new variable v = jy, where j = 

The contribution from Matsubara frequencies with 1, 2, 3, .... Thus, we obtain 
/ > 1, ^n~^Va,T), is more cumbersome. We will find 

its low-temperature asymptotic behavior perturbatively. ^ _ 21k_2 y^n _|_p2(^2N-3/2 £ / dyg-^ ( Ci — — ) 

With this purpose we expand Zd in Eq. ([52|l in powers of 2LUp ^-^ ' ^-—^^ j J \ J 

a small parameter xi ^gj^-j 

„A f „^ ^-2^ After the integration with respect to v, this leads to 

where the impedance of infrared optics, Z^, is defined in 

(62) 

Eq ([35]) and expanding in powers of pQ = pul and performing the 

summation in I we obtain 

Now, we substitute Eq. ([55)1 in Eq. ([M]) with omitted zero- 
frequency contribution and expand J-'^^^^ (a, T) in powers g 
of keeping only the first order term: -f — p^K^ 



2' 



1 ^ e'J"^ ^£2^(1 + 6-^") 



T^^^'\a,T)=Tj:'^'\a,T) + ^Y.^iii (58) +0(pV)}. (63) 



y dye' 



-y + (iZ^iiQ^c) , Qi - yZi{iC,iQc) 



In the asymptotic limit k — > 0, Eq. (|63p results in 



The substitution of this expression into Eq. ([591) leads to 



^TM ^TE J 7(r) luK , , 7(T) /I", 

- ^«IL^_,;(3) + IL^0 - . (64) 

-fO(a::;), ujp n Up \k/ 

where 

^TM = [y + CiZiiiQ^c)] 

X [(e-^ - 1) {y' + Cfzn^CA)) T^^'\a,T) = :F^'^'\a,T) + M^iZl i^i^^O) 

-2{e-y + l)0yZ,it0n,)], 27raa;p 

Ate = [0 + yZSCifi.)] + O [Z,itOn,)i,,if] + O f . (65) 

X [{e-y 1) id + y'Zfi^Q^,)) 2™ 

—2 {e.~y + 1) C,iyZi{iC,iQc)\ ■ ^ow we are in a position to find the asymptotic ex- 
pression at small k for the free energy computed using 
Here Tf-'^\a,T) is the free energy of the fiuctuating the impedance of the Drude model. First we add the zero- 
field, computed by using the impedance of infrared optics frequency term !F^^~^\a, T), defined in Eq. (|55p . to both 



10 



V. B. Bezerra et al.: Thermal correction to the Casimir force and radiative heat transfer 



sides of Eq. ([65| . On the left-hand side of this equation 
the quantity J-oio-, T) is obtained. On the right-hand side 
of Eq. (j65[) we add and subtract the zero-frequency term 
of the free energy computed using the impedance of the 
plasma model. This term is obtained from Eqs. (|36p - (|55|) : 



In 



1 + py 



ydy I In (l - e ^) 



(66) 



The term (|66|) . together with the first term on the right- 
hand side of Eq. gives us the free energy J^i{a,T) 
computed using the impedance of infrared optics. As a 
result, from Eq. ((65|) it follows: 

TD{a,T) = T^{a,T) + ATf=''^ (a, T) 
+ — C(3) + O [ZSQO^)iuil] 



kBTj{T) ^fl 



(67) 



where 



AtI^-°'> {a, T) = (a, T) - (a, T) (68) 



j ydy I In (1 - e ^) 




Taking into account that k ^ T and at low temperatures 
7(r) - T^, Z^iiCiHc) ^T,ii'^"f ^ T^, we arrive at the 
conclusion that not only the last three terms on the right- 
hand side of Eq. (j67p vanish when temperature vanishes, 
but also their derivatives with respect to temperature van- 
ish. 

Using Eq. (|50p . we can find the asymptotic behavior 
of the free energy and entropy of the fluctuating field at 
low temperatures in the case that the metal is described 
by the impedance of the Drude model. Keeping only the 
main terms of order T and In T in Eq. ([67|) [recall that 
according to Eq. ([49]) Ti{a, T) - E{a) ~ T^], we obtain 



TD{a,T) = E{a) + ATf-''\a,T) 
, fcBT7(r) In^^^g^ 



Soia^T) 



dAT, 



(69) 



p 

(i=0) 



(a,r) 



fcBC(3) l{T) Teff 




167ra2 I ./o 
27r2a2 



Expanding the integrand on the right-hand side of the 
second equality in Eq. (j69p in powers of p and integrating 
with respect to y, we arrive at 



SD{a,T) 



fcsC(3) 



27ra2 
fcBC(3) 7(J^) J^ff 



p[l-6p + 0(p2)] 



(70) 





1" 


In f 27r— I 






+ 2 



To carry out the thermodynamic test of the impedance 
of the Drude model, we consider T approaching zero and 
find the following value of the entropy: 



SD{.a,Q) 



16na 



C(3) 



ydy In 
ksm 



py 



py 



p[l-6/? + 0(p2)] >0. (71) 



dT 



As is seen from Eq. ([7T|) , the entropy of the fluctuating 
field at zero temperature takes a nonzero positive value. 
This value depends on the parameters of the system, i.e., 
on the separation distance a, and on the plasma frequency 
ujp. The latter participates through the definition of p (re- 
call that for metals without impurities described by 
the Casimir entropy at T = is negative [H^HO] ). What 
this means is that we have a violation of the third law of 
thermodynamics, the Nernst heat theorem [7^ . 

In the above calculations we have considered the en- 
tropy associated with the fluctuating field. In other words, 
only the distance dependent part of the total free energy 
was considered. However, inclusion of the self-energies can- 
not invalidate our conclusion on the violation of the Nernst 
heat theorem. The reason is that the entropy of a fluctu- 
ating field in Eq. (|7ip depends on the separation distance, 
whereas the entropies due to self-energies are separation 
independent. As a consequence, the entropy of a fluctuat- 
ing field and the entropies of matter fields cannot cancel 
each other and must satisfy the Nernst heat theorem sep- 
arately. 

Thus, the Leontovich impedance of the Drude model 
is thermodynamically inconsistent and cannot be used 
in combination with the Lifshitz formulas, Eq. ([T]) and 
Eq. p2|) . to calculate the thermal Casimir force. In the 
next section it is also shown that the thermal Casimir 
force computed using the impedance of the Drude model 
is in disagreement with the experimental data of |17j . 
Thus, there are both theoretical and experimental evi- 
dences against the use of this impedance function in the 
theory of the thermal Casimir force. 
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6 Alternative results for the thermal 
correction to the Casimir pressure computed 
with different impedance functions 

In [571158] the contribution from the transverse electric 
electromagnetic waves to the thermal correction 
AP{a, T) was computed using the real frequency axis for- 
malism and the impedance function of the normal skin 
effect Zjs[{u}) defined in Eq. ((32)) . The computation was 
performed for Au plates at room temperature with static 
conductivity ctq = 3 x 10^^ s~^ and relaxation time r = 
1.88 X 10""'^"' s. As was noted in Sec. 3, the impedance func- 
tion p2|) can be obtained from the dielectric function in 
the region of the normal skin effect, ejv(tj), using Eq. (|26p 
where the static conductivity is connected with the plasma 
frequency ujp by Eq. (134]) leading to ujp w 9.3 eV. In [571158] 
the relatively large contribution to the thermal correction 
from the TE waves for plates separation a = 1 /im was ob- 
tained using the impedance function ([5^ . It is about 30 
times larger than the corresponding correction for ideal 
metals. This was explained by the increased role of the 
TE EW at low frequencies. According to [SS] the predicted 
increase of the thermal correction at a = 1 /xm is in con- 
flict with previous theoretical work |;37 that predicts even 
several times larger magnitudes for the thermal correction 
calculated using the dielectric permittivity e ^ cj^^ at low 
frequencies. Paper |59j concludes also that the prediction 
of [57,58^ is consistent with the experimental data of [9\ 
at a = 1 /Ltm whereas the prediction of 37^ is excluded by 
that experiment. 

Here we repeat the numerical computations of [571158] 
for the contribution of the TE mode to the thermal cor- 
rection to the Casimir force, using the impedance ([5^ . 
We also compute the contribution from the TM mode 
and study the role of EW and PW in both contributions. 
The computations are performed at different separations 
using both formalisms along the imaginary and real fre- 
quency axis (presented in Sees. 2 and 3) with practically 
coinciding results. Although the thermal correction pre- 
dicted from the impedance ([32]) is consistent with the long- 
separation experiment [9] at 1 ^m, it is shown to be ex- 
cluded by the measurement of the Casimir force at shorter 
separations by means of the micromechanical torsional 
oscillator [17,18,19.20] . Because of this, we also discuss 
several other forms of impedance function and find those 
consistent with experiment. 

For the purpose of comparison with experiment (which 
is in fact consistent with the theoretical results for Po{a) 
at T = [T7]), here we use in computations slightly dif- 
ferent value of the plasma frequency ujp = 9.0 eV. The 
thermal correction in the framework of the imaginary fre- 
quency axis formalism was computed by Eqs. ([2]), ([5]) 
and ([9]) for two similar Au plates at T = 300 K. The re- 
sults are presented in Table 1 where column 1 contains 
the values of the separation distance between the plates 
and column 2 the values of the thermal correction. Pre- 
cisely the same values were computed using Eq. ^TE\\ of 
the real frequency axis formalism. Note that although 
the magnitude of the thermal correction increases with 



the decrease of separation, the relative thermal correc- 
tion becomes smaller at shorter separations. The same 
holds for metals described by the dielectric permittivity 
of the plasma model [34,35J. In column 3 of Table 1 the 
ratios of the thermal correction from column 2 to those for 
ideal metal (denoted by Z\P^e) presented. Columns 4 
and 5 contain the relative contributions to the thermal 
correction from the TE EW and PW, respectively, com- 
puted by using Eqs. ([20]) and ([2T]) . In columns 6 and 7 
the relative contributions from the TM EW and PW, re- 
spectively, are presented. Here and below we perform all 
computations at a > 200 nm in order to remain well in- 
side the application region of the impedance approach. 
For example, at a = 200 nm the characteristic frequency 
of the Casimir force is ujc = cj (2a) « 7.5 x 10^^ rad/s and 
Zn{^c) ~ 1.4 X 10^^ ^ 1. As is seen from columns 2 and 3 
in Table 1, the magnitudes of the total thermal correction 
computed with the impedance function ()32|) are rather 
large. At separation distance of 1 the thermal correc- 
tion for Au plates is found to be almost 16 times larger 
than for the plates made of ideal metal. This ratio quickly 
increases with the decrease of separation. At a separation 
of 200 nm it is as large as 5900. By the summation of the 
values presented in columns 4, 5 from one hand and 6, 7 
from another hand, one finds that the dominant contribu- 
tions to AP are given by the TE mode, whereas the rela- 
tive contributions from the TM mode are negligibly small. 
The largest value of the latter achieved at a = 1 /im is 
equal to 0.05. Comparing columns 4 and 5 from one hand 
and 6 and 7 from another, we can conclude that the dom- 
inant contribution to Z\Pte and APtm is given by the 
EW. If one considers only the contribution from the TE 
mode discussed in [ST], one obtains Z\Pte/^^te ~ ^9.6 
and 1.2 x 10^ at separations a = 1 /im and 0.2 /im, respec- 
tively. 

Now we compare the obtained magnitudes of the ther- 
mal correction in column 2 of Table 1 with the experiment 
117] . For this purpose we need to obtain the total magni- 
tudes of the Casimir pressure at temperature T = 300 K. 
The magnitudes of Po(a) can be most simply computed 
using Eqs. ([9]) and ([27]) . However, use of the impedance 
function ([32]) leads to incorrect values of Pq (a) . As an ex- 
ample, using ([32]) one obtains \Po\ = 716.5 and 144.6 mPa 
at separations 200 nm and 300 nm, respectively. At the 
same time, the conventional magnitudes of Pq at these 
separations, obtained by different authors [76l[77l[78] are 
507.5 and 113.6 mPa, respectively, in drastic contradiction 
with the above values obtained by the use of impedance 
([32]) . The reason is that the dominant contribution to Pq 
is given by the frequency region around the characteris- 
tic frequency coc which belongs not to the region of the 
normal skin effect, where the impedance function (|32p is 
appropriate, but to the region of infrared optics. On the 
contrary, the magnitude of the thermal correction (|18p 
with the impedance function (|32|) is determined by much 
lower frequencies where ([32]) is applicable. 

To compare the thermal correction in column 2 of Ta- 
ble 1 with experiment, we add it to the conventional mag- 
nitudes of the Casimir pressure at T = specified above. 
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As a result, at separations 200 and 300 nm one obtains the 
magnitudes of the thermal Casimir pressure equal to 519.6 
and 116.4mPa, respectively. These should be compared 
with respective measured magnitudes of the Casimir pres- 
sure at the same separations equal to 508.1 and 114.7 mPa. 
The differences between the above theoretical and exper- 
imental values are 11.5 and 1.7niPa. They lie outside the 
boundary of the half- width confidence interval determined 
at 95% confidence (at separations 200 and 300 nm the lat- 
ter is equal to 8.6 and 1.6 mPa, respectively). Thus, the 
impedance function (|32p is not consistent with the mea- 
surement of the Casimir pressure at short separations by 
means of micromechanical torsional oscillator [TT l lTS l fTQ l 
[20]. 

The impedance Zn is applicable only at low frequen- 
cies specific for the normal skin effect. As was discussed in 
Sec. 3, the impedance of the Drude model ((33)) provides 
smooth interpolation between the regions of the normal 
skin effect and infrared optics. At all frequencies uj ujp 
the impedance p3|) can be represented in the form (|32p if 
one replaces ctq on the right-hand side of ((32|) by the ac 
conductivity defined as 

a{u;) = (72) 

1 — ITLO 

It is easily seen that the impedance of the Drude model 
([55]) leads to even larger magnitudes for the thermal cor- 
rection to the Casimir pressure at short separations, than 
the impedance (|32|) . At the separation a = 0.2 /^m the 
magnitude of the thermal correction computed with the 
impedance ([55]) is equal to — 14.7 mPa and its ratio 
to the thermal correction for ideal metals is 7200. At a = 
1/xm one obtains |Z\P| = 0.0402 mPa and Z\P/Z\F™ = 
20. If one considers only the contribution from the TE 
mode, one obtains Z\Pte/^Jte = 37.4 (close to 36.5 
obtained in [58j for a bit different value of uip — 9.3 eV 
whereas in our computations here we use LUp — 9.0 eV). 
However, at a w 200 nm this ratio achieves the value 
Z\Pte/^-P|e ~ 14.3 X 10^. As in the case of impedance 
function ((3^ , the increase of the thermal correction is due 
to the contribution from the TE EW. 

In the same way, as before, it can be shown that the 
large thermal correction predicted by the impedance p3|) 
is excluded experimentally. Thus, the theoretical magni- 
tudes of the Casimir pressure obtained by adding the ther- 
mal correction to Pq are equal to 522.15 and 116.99mPa 
at separations 200 and 300 nm, respectively. By compar- 
ing this with experiment we get the respective differences 
14.06 and 2.32 mPa which are far outside the boundary of 
the 95% confidence interval presented above. [Notice that 
the use of the Drude dielectric function from Eq. p3|) with 
reflection coefficients ([22]) to compute the thermal Casimir 
force is inconsistent not only with experiment [9J but with 
experiments [l6l[T7p8l[T9pO] as well.] 

Bearing in mind that in the separation range from 
200 nm to 1/im the characteristic frequency i7c ~ c/(2a) 
belongs to the region of infrared optics, it is reasonable to 
calculate the thermal correction to the Casimir pressure 
using the impedance Zi defined in Eq. (|3ip . Note that the 
dielectric permittivity ei(w) with the reflection coefficients 



([22]) was first used to calculate the thermal Casimir pres- 
sure in [34j[35] . It was shown that the dielectric permittiv- 
ity of the plasma model leads to small thermal corrections 
to the Casimir pressure at short separations in qualitative 
agreement with the case of ideal metals. At large sepa- 
rations it leads to the same result as for ideal metals in 
accordance to the classical limit [73)1 74 j. 

The disadvantage of the impedance function (pij) is 
that it has zero real part at nonzero frequencies, i.e., does 
not take relaxation processes into account. At the same 
time the tabulated optical data for the complex index of 
refraction n in the region of infrared optics [75] lead to a 
nonzero real part of the impedance function which, how- 
ever, is much less than the magnitude of the imaginary 
part. Unfortunately, the optical data at low frequencies 
are not available. In addition, these data are burdened by 
large errors and uncertainties. Because of this, any the- 
oretical result for the functional form of the impedance 
function is of much value. Preserving only the first ex- 
pansion orders for both real and imaginary parts, a more 
accurate complex impedance in the region of the infrared 
optics can be approximately presented in the form [5T] 

Zpico) = Cuj^ - i (73) 

where the constant C = 0.004 eV~^ was chosen to provide 
the best mean fit to the optical data for the real part of the 
impedance within the frequency region from 0.125eV to 
5 eV. The impedance ((73)) disregards interband transitions 
but the real part of it takes into account electron-electron 
collisions which occur rather seldom in the region of in- 
frared optics. 

In Fig. 1 the imaginary (a) and real (b) parts of the 
impedance ()73p are shown by the sohd lines as functions 
of the frequency. In the same figure the imaginary (a) and 
real (b) parts of the impedance are plotted by dots using 
the expression Zp{uj) — \/n{Lj') and the tabulated optical 
data for n{uj) 75J. As is seen in Fig. lb, in the region 
w < 1.5 eV the analytic expression for the imaginary part 
of the impedance ([73]) is in rather good agreement with 
data. This expression gives the major contribution to the 
impedance along the imaginary frequency axis 

ZM) = -Ce' + (74) 

Then one can conclude that the approximation ()74p is well 
adapted for the computation of the thermal correction to 
the Casimir pressure. 

The computations were performed for Au plates by us- 
ing both the imaginary and real frequency axis formalisms 
discussed in Sec. 2 with practically coinciding results for 
the total thermal correction and also for separate contri- 
butions from the TE and TM modes. When using the real 
frequency formalism, some care is necessary in the numer- 
ical integration of Eqs. ([20)) and ()2ip . because, due to the 
smallness of the real part of Zi, there exist narrow reso- 
nances in the spectrum of the thermal correction, corre- 
sponding to the eigenfrequencies of a dissipationless cavity 
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Fig. 1. Imaginary (a) and real (b) parts of the impedance of 
infrared optics as functions of frequency. Solid lines are for the 
analytical expression (73), dots correspond to the tabulated 
optical data. 



with impedance Zi{uj). These considerations apply espe- 
cially to the computation Z\Ptm,eWj because its spectrum 
has resonances also at thermal frequencies. These results 
are presented in Table 2. In column 2 the magnitudes of 
the total thermal correction to the Casimir pressure at dif- 
ferent separations are listed. Column 3 contains the ratios 
of the total thermal corrections for Au plates to those for 
ideal metals. In column 4 the relative contributions from 
the TE EW to the total thermal correction are shown. Col- 
umn 5 contains the relative contributions from the TE PW 
to the total thermal correction. Columns 6 and 7 contain 
the relative contributions to the total thermal correction 
from the TM EW and TM PW, respectively. As is seen 
from columns 2 and 3 in Table 2, the impedance of infrared 
optics leads to much smaller thermal corrections than the 
impedance of the normal skin effect in qualitative agree- 
ment with the case of ideal metals. From columns 4-7 it 
follows that the TE and TM modes lead to qualitatively 
similar contributions. The role of EW is also not so pro- 
nounced as it was for the impedance of the normal skin 
effect. 

The total thermal Casimir pressures obtained by the 
summation of the thermal corrections in column 2 of Ta- 
ble 2 and zero-temperature pressures are consistent with 
all available experimental data. By way of example, at 
a = 200 and 300 nm (where for the impedance of the nor- 



mal skin effect the theoretical results were excluded at 
95% confidence) the differences between theory and ex- 
periment ^17,jl8| are now equal to -0.61 and -1.1 mPa, 
respectively, i.e., well inside the 95% confidence interval. 
Note that the thermal Casimir force computed using the 
dielectric permittivity pT|) is also consistent with data 
\2-7\. The addition of a small imaginary part to Si{uj), aris- 
ing from real part of the impedance in Eq. (|73p . leads 
to only minor changes in the computation results which 
remain consistent with experiment. 

Thus, the impedance of infrared optics (|74p leads to 
reasonable results when applied to the thermal Casimir 
pressure. In the following sections we are going to apply 
the impedance method to the computation of the radiative 
heat transfer across a vacuum gap between two parallel 
surfaces at different temperatures. As is shown in Sec. 7 
below, this process is mostly determined by the real part 
of the impedance function which is modelled not enough 
precisely in Eq. (|73p (see Fig. lb). Bearing in mind the 
computations of the radiative heat transfer, we introduce 
one more model impedance of infrared optics by the equa- 
tion 




uj<(3, 

(3<uj< 0.125 eV, 
w > 0.125 eV 



(75) 



Here Y{u!) stands for the tabulated optical data rep- 
resenting the real part of the impedance [75j which are 
available at a; > 0.125 eV. The value of the constant B = 
0.00389 is chosen such as to have a smooth transition be- 
tween the optical data and their extrapolation to lower 
frequencies. The unknown parameter /3 fixes the value of 
frequency where the behavior of ReZj, as given by the 
optical data, is smoothly connected with the asymptotic 
behavior at low frequencies. The upper bound of P is de- 
termined by the fact that the optical data are available at 
uj > 0.125 eV. In the frequency region of infrared optics 
the real part must be much smaller than the imaginary. 
The reason is that at these frequencies electrons are almost 
free and electric current is pure imaginary. The small real 
part of the impedance describes minor distortions in the 
vibrational motion of electrons. Then it is reasonable to 
impose the lower constraint on /3 as follows: 0.08 eV < /3. 
This constraint is very conservative because it allows the 
real part of the impedance to become as large as one half 
of the imaginary part (a larger real part of the impedance 
in the region of infrared optics is evidently inadmissible). 
The real part of the impedance ([75]) as a function of fre- 
quency is shown in Fig. 2, where line 1 corresponds to 
/3 = 0.125 eV and line 2 to /3 = 0.08 eV. The region be- 
tween lines 1 and 2 is allowed. 



7 Heat transfer across an empty gap 

We consider now the case of two semi-infinite plane par- 
allel metallic plates at different temperatures Ti > T2, 
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Fig. 2. Real part of the impedance of infrared optics given by 
the tabulated optical data (dots) versus frequency with differ- 
ent extrapolations to low frequencies (the region between the 
solid lines 1 and 2). See text for further discussion. 



separated by an empty gap of width a. We are interested 
in estimating the heat transfer between the plates. This 
problem was first studied long ago by Rytov [79] , and later 
on it was reconsidered by Polder and Van Hove [5D], by 
Loomis and Maris [61] and by Volokitin and Persson [Sj • 
Recently, it was studied by one of us ^3] in the framework 
of the surface impedance. The approach followed by Polder 
and Van Hove is closely related to Lifshitz theory of the 
van der Waals interactions between macroscopic bodies, in 
that heat transfer is regarded as occurring via fluctuating 
electromagnetic fields radiated by the two plates, whose 
sources are the random thermal electric currents that are 
present inside the plates. The fluctuation-dissipation the- 
orem is then used to determine the statistical properties 
of the currents, from which the correlation properties of 
the emitted electromagnetic fields are subsequently de- 
rived. In the derivation, Polder and Van Hove made the 
following assumptions, that are also at the basis of Lifshitz 
theory: i) the wavelengths of the electromagnetic fields in- 
volved should be large compared to atomic distances, so 
that the fluctuating flelds can be well described by means 
of the classical macroscopic Maxwell's equations; ii) the 
electric currents at distinct points inside the plates are un- 
correlated; iii) the media are isotropic and nonmagnetic, 
and such that their electromagnetic properties can be de- 
scribed by means of a complex dielectric permittivity e{uj) 
that depends only on the frequency; iv) the system is sta- 
tionary in time, and each plate is in local thermal equi- 
librium. The resulting expression for the power (per unit 
area) S of heat transfer from plate one to plate two, was 
found using the average value of the Poynting vector in 
the gap between two plates [61]: 



S 



Ah 
1^ 



1 



dk \ k \ k 



2\2ik^a\ 



(76) 



^exp{huj/kBTi) — 1 exp{hLo/kBT2) — 1 
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X 



TM 



where 



XTE = |(fc. + S«)(fc, + S^^J) 



,(2)1 



=(1) 



(2)-jg2ifc,a|2 



XTM = |(eW A:,-f.«)(£(2) A:, + .(2)) 

_(e(i) A:,-s«)(£(2) fc,-s(2))e2»fc.a|2^ 

and 

s("Ht^,fc_L) = ■y/e(")(w)t^7c2-fc2 =i/c(")(cj,fc_L) (77) 

(note that our variable k± is denoted by q in [5^61j ). We 
now decompose S as the sum of the contributions from 
PW and EW: 

S = S'pw + 'S'ew • (78) 

It is not hard to verify, starting from Eq. (j76p . that 5pw 
and Sew can be expressed in terms of the dielectric re- 
flection coefficients (l24l) as 1321: 



S- 
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It is to be noted that, for r, 

,(2 



(2) 



= 0, Sew — Oi while for 
ra' =0, and T2 = the expression for S'pw reduces to the 
well known Kirchhoff 's formula for the flux of radiation 
from a surface with reflection coefficients = ra\ at 
temperature T = Ti: 
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47r2 c2 



duj ■ 



exp{huj/kBT) — 1 



dp p 



a=TE,TM 



(1- 



(81) 



where p — k^ c/uj. 

In [33] the power of heat transfer per unit area S was 
computed within a general theory of electromagnetic fluc- 
tuations for metallic surfaces, based on the concept of sur- 
face impedance. This approach is closer to that originally 
followed by Rytov [79 , in that the starting point of the 
theory is an expression for the correlators of the electric 
and magnetic fields. However, according to the dictate of 
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Fig. 3. Plots of radiative heat transfer between two Au plates 
at temperatures Ti — 320 K and T2 = 300 K, as a func- 
tion of separation, according to Lifshitz theory for e = ed 
(short-dashed line), and to impedance theory for three differ- 
ent choices of impedance: Z = Zn (long-dashed line), Z = Zd 
(dotted line) and Z — Zt (band between the solid lines 1 and 
2). See text for further explanation. 



impedance theory, in the correlators are given only 
outside the metal, while no consideration is made of ei- 
ther the fields, or the electric currents in the interior of 
the metal. The resulting expression for the heat transfer 
that was found in 1331 is: 



dio 
1 



ReZ(i)(cj)Re^(2)(^-) 
1 



exp{hu}/kBTi 



dk I k \ \k. 



- 1 

2ik, 



exp{h,u)/kBT2) - 1 
1 1 



,(82) 



/o \Bte Btm , 

where the quantities Bj-^/tm are defined as: 

Bte = \{l + k, c Z'-^^ /uj){l + k, c Z'-^^ /lu) (83) 

- (1 - k,cZ^^'^/u;){l - k,cZ^^^/u;)exp{2ik, a)p , 
Btm = |(fc. c/c^ + Z(i))(A:, c/lu + Z^^)) (84) 

- {k^ c/uj - Z^^^){k, c/lu - exp(2i/fc^ a)p . 



If we separate the PW and EW contributions to Eq. 
it is easy to verify that Eqs. (|75 |l - ([5D)) are equivalent to 
Eqs. in the impedance theory of heat transfer, 

provided that the reflection coefficients are taken to be 
those of impedance theory in Eq. ([28]) . 



8 Numerical results for heat transfer and 
emittivity 



impedance function, at angular frequencies lu of the or- 
der of kBT/h. It is also important to observe that at 
small separations a, S receives a large contribution from 
the EW. Since, as seen from Table 1, thermally exited 
EW are of great importance for the determination of the 
thermal correction AP{a, T) to the Casimir pressure, it 
is clear that a measurement of S would provide an in- 
dependent verification of the validity of the impedance 
function, used in the evaluation of Z\P(a, T). We have es- 
timated numerically the heat transfer S for three choices 
of the impedance functions, Z^, Zjq and Zt, that were 
discussed in Sec. 6. In Fig. 3, we show plots of the ra- 
diated power for two plates of Au as a function of the 
separation a for Ti = 320 K and T2 = 300 K. The three 
lines are for the standard Lifshitz theory with the Drude 
dielectric function in Eq. ([55]) (short-dashed line), for the 
impedance theory with the impedance Z^ of the normal- 
skin effect in Eq. (l32|) (long-dashed line), and again for 
the impedance theory, but this time with the impedance 
Zd corresponding to the Drude model in Eq. p3|) (dotted 
line). The band between the solid lines 1 and 2 is related 
to the impedance function Zt in Eq. ([75l) . with the upper 
and lower boundaries corresponding to (3 — 0.08 eV and 
P = 0.125 eV, respectively The hne for Z = Zp in Eq. ^ 
is not displayed, because it reproduces ReZp inaccurately 
and the corresponding values for S are over two orders of 
magnitude smaller than those for, say, Zt. Note that the 
dielectric permittivity approach with e corresponding to 
Zt leads to almost the same results as are presented by 
the band between the solid lines 1 and 2 in Fig. 3. The 
only difference is that at short separations the lines 1 and 
2 are a bit shifted towards smaller values of S, while pre- 
serving the same asymptotic values at large separations. 
As we see, Lifshitz theory using the dielectric permittiv- 
ity Ed, as well as the impedance theory for Z = Zo and 
Z = Zm, both lead to values for S that are several times 
larger than those implied by the impedance Zt or respec- 
tive dielectric permittivity, for separations around or less 
than half micrometer. To a large extent, these large differ- 
ences are due to considerably different contributions from 
the TE EW in the various models. In Fig. 4, we show a 
plot of the relative contribution to S from TE EW, for 
the models considered in Fig. 3. As we see, the TE EW 
play an important role in the entire range of separations 
considered. 

Besides heat transfer, another interesting quantity to 
consider is the total emittivity e{T) of the metal, defined 
as: 

e(T) = . (85) 

Here, ^(T) is the total flux of radiation from a unit surface 
of a metal defined in Eq. ((8T|) . while '1'bb{T) is the flux 
from a black body. According to Stefan law, the latter 
c^uantity is equal to 



(86) 



It should be noted that, according to Eq. ([5^ . the value 
of S is very sensitive to the real part ReZ(a;) of the 



where = 5.6704 x 10"** W m^^K""*. For a polished sur- 
face of gold at T = 295 K, the tabulated value is e w 0.02. 
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Fig. 4. Plots of the relative contributions of TE EW to the 
total radiative heat transfer between two An plates at temper- 
atures Ti — 320 K and T2 — 300 K as a function of separation. 
Lines are notated as in Fig. 3. 

In Table 3 we report the calculated values for e, corre- 
sponding to the same models as considered in Fig. 3. The 
interval of values for e, appearing in the fifth column in the 
case oi Z = Zt, is related to the values of /3 considered in 
Fig. 3, with the lower and upper values of e corresponding 
to (3 — 0.125 eV and /3 = 0.08 eV, respectively. As is seen 
in Table III, the measured emittivity is better described 
by the impedance of the normal skin effect and by the 
generalized impedance of the infrared optics Zt. However, 
any definite conclusion is impossible without information 
on the precision of emittivity measurements. 

9 Conclusions and discussion 

In the first part of this paper we have implemented the 
thermodynamic test for two different choices for the sur- 
face impedance of a metal, that are used in the theory of 
the thermal Casimir force. By making analytic perturba- 
tion expansions in powers of small parameters, we have 
obtained the asymptotic expressions for the free energy 
and entropy of a fluctuating field at low temperatures. 
This was done by using the Leontovich impedances of the 
infrared optics Zi and of the Drude model Zd . The Leon- 
tovich impedance of the infrared optics, extrapolated to 
low frequencies, withstood the thermodynamic test. The 
obtained asymptotic expressions for the free energy and 
entropy of a fluctuating field are found to be thermody- 
namically consistent. In particular, entropy becomes zero 
when temperature vanishes, i.e., the Nernst heat theorem 
is satisfied. On the other hand, the Leontovich impedance 
of the Drude model was shown to be thermodynamically 
inconsistent. In the case of metals with perfect crystal 
lattices the entropy at zero temperature was found to be 
positive and depending on the parameters of the system 
in violation of the Nernst heat theorem. 

The above conclusions provoke two questions on how 
to correctly apply thermal quantum field theory in Mat- 



subara formulation to real materials. Both the result of 
this paper and of |39p40j on the thermodynamic incon- 
sistency of Zu and Ed, respectively, in the theory of the 
thermal Casimir force were obtained by using the ideal- 
ization of perfect crystal lattice of a metal. In the pres- 
ence of impurities there is a nonzero residual relaxation 
at T = 0, i.e., 7(0) 7^ 0. As a result, at very low temper- 
atures the first Matsubara frequencies may become less 
than 7(0) and the entropy jumps steeply to zero. For met- 
als with impurities, described by the dielectric permittiv- 
ity of the Drude model, the vanishing of entropy at T = 
was demonstrated in [4T|l42j. This, however, does not solve 
the problem arising for metals with perfect crystal lattices. 
Such metals have a nonzero relaxation 7(r) at nonzero T 
and are commonly used as the basic model in the theory 
of electron-phonon interactions. For perfect crystal lattices 
with no impurities the Nernst heat theorem is proved in 
the framework of quantum statistical physics |72j , and the 
violation of this theorem by the entropy of a fluctuating 
field is a problem of great concern. It is our opinion that 
the violation of the third law of thermodynamics by the 
Casimir entropy calculated using ed and Zjy warns about 
the inapplicability of the Drude model in the theory of the 
thermal Casimir force. 

Another question to discuss is the physical meaning 
of the zero Matsubara frequency in the Lifshitz formula, 
Eq. ([1]), or of low frequencies in the equivalent form of it, 
Eq. (jl2p . expressed in terms of real frequenceis. Should the 
analytic expression for the impedance of a metal which is 
valid for frequencies around the characteristic frequency 
f2c, be extrapolated without modifications to quasistatic 
frequencies, as we did in Sec. 4, where we dealt with the 
impedance of infrared optics? Or, alternatively, should 
one use different impedance functions within different fre- 
quency regions in accordance with their applicability con- 
ditions? The latter approach was in fact used in Sec. 5 
because the impedance of the Drude model coincides with 
the impedance of the infrared optics in the region of in- 
frared frequencies around f2c and with the impedance of 
the normal skin effect in the region of quasistatic frequen- 
cies. Our results demonstrate that in spite of being rather 
natural to use different impedance functions in accordance 
with the frequency regions of their applicability, and not 
use any extrapolation, the actual situation is not so simple. 
As is shown in Sec. 4, the extrapolation of the impedance 
of infrared optics to zero Matsubara frequency satisfies the 
thermodynamic test, whereas the use of the Drude model 
impedance in Sec. 5, coinciding with the impedances of the 
normal skin effect and infrared optics in the appropriate 
frequency regions, violates thermodynamics. This should 
be compared with the results of [44ll45p 46 p 47 | l48] devoted 
to the Casimir interaction between two dielectrics and be- 
tween metal and dielectric. In both cases the account of an 
actual dielectric response at quasistatic frequencies (i.e., 
the account of nonzero dc conductivity) results in contra- 
diction with thermodynamics, whereas the extrapolation 
of the dielectric behavior at high frequencies to zero fre- 
quency satisfies the thermodynamic test. This leads us to 
argue that the response function of both a metal and a 
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dielectric to a real external electromagnetic field of very 
low, quasistatic, frequency is not related to the physi- 
cal phenomenon of dispersion forces determined by the 
electromagnetic fluctuations of high frequencies. In terms 
of physical processes occuring at different frequencies the 
above discussed extrapolation would imply that in metal 
bodies the fluctuating electromagnetic field creates only 
pure imaginary currents related to the frequency region 
of infrared optics. However, if the characteristic frequency 
belongs to the region of infrared optics, real currents with 
typical frequencies of the normal skin effect cannot be cre- 
ated by the fluctuating field. The deeper understanding of 
these gueses may go beyond the scope of the Lifshitz the- 
ory. 

In the second part of this paper we have performed a 
comparative phenomenological investigation of the ther- 
mal Casimir force between two parallel metal plates and 
of the radiative heat transfer which occurs between such 
plates when they are kept at different temperatures. Both 
phenomena are of the same physical nature because they 
are caused by electromagnetic fluctuations. As was dis- 
cussed above, the problem of the thermal Casimir force 
meets with difficulties, and different controversial appro- 
aches to its resolution were proposed in the literature. 
These approaches are based on the use of various dielec- 
tric functions (the dielectric permittivities of the Drude 
and of the plasma model) or, alternatively, different forms 
of the Leontovich surface impedance. The selection be- 
tween the approaches is done by the comparison of the 
obtained results with the requirements of thermodynam- 
ics and with the experimental data. In this paper we have 
demonstrated that the use of the impedance functions 
Z]\r{ui) and Zd{uj) (constructed using the dielectric per- 
mittivities of the normal skin effect and of the Drude 
model) to calculate the thermal correction to the Casimir 
force leads to contradiction with experiment at separa- 
tions of a few hundred nanometers. Recall that earlier the 
use of the dielectric permittivity of the Drude model to 
calculate the thermal correction was also shown to be in- 
consistent with experiment. At the same time, both the 
dielectric permittivity of the plasma model and the cor- 
responding impedances related to the region of infrared 
optics are consistent with experiment. 

The radiative heat transfer between two plates was 
previously studied using the dielectric permittivity of the 
Drude model [S^USDllBlj and the related impedance Zu{(^) 
[33] . These approaches, however, lead to a contradiction 
with experiment in the case of the thermal Casimir force. 
Because of this, it is of much interest to investigate the 
radiative heat transfer using the impedance of infrared 
optics. Lack of reliable experimental data for the power 
of heat transfer and resulting uncertainty with the experi- 
mental confirmation of the computations in [32 l [33 l [60 j l6T] 
add importance to this aim. 

Bearing in mind that the radiative heat transfer is very 
sensitive to the real part of the impedance function we 
have constructed the new impedance Zt in the region of 
infrared optics. The real part of this impedance takes into 
account the tabulated optical data for the complex index 



of refraction extrapolated to low frequencies in accordance 
with general theoretical requirements. The power of heat 
transfer calculated with this impedance is several times 
less than previous predictions at separations of a few hun- 
dred nanometers. These large differences are mainly ex- 
plained by different contributions from the TE EW in the 
various models of a metal. 

Both physical phenomena of the thermal Casimir force 
and of the radiative heat transfer are finding prospec- 
tive applications in nanotechnology. This makes urgent to 
carry out new precise experiments in order to find what 
characterization of real metals is most adequate for the 
description of electromagnetic fluctuations. 
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Table 1. Thermal correction to the pressure between Au 
plates at r = 300 K (column 2) and different contributions to 
it as a function of separation computed using the impedance 
of the normal skin effect Zn- See text for further discussion. 



a (/im) 


AP (mPa) 


AP 


^-Pte,ew 


^-Pte,pw 


^-Ptm,ew 


^-Ptm,pw 




AP 


AP 


AP 


AP 


0.2 


-12.1 


5.9 xlO-^ 


0.998 


-9 xlO"^ 


2 xlO"^ 


-2 xlO"* 


0.25 


-5.4 


2.6 xlO^ 


0.997 


7 xlO"^ 


3.6 xlO"^ 


-2 xlO"* 


0.3 


-2.8 


1.4 xlO^ 


0.995 


-3 xlO"*^ 


5 xlO"^ 


-3 xlO""* 


0.35 


-1.6 


7.8 xlO^ 


0.994 


1.5 xlO"'' 


6 xlO"^ 


-3 xlO"* 


0.4 


-0.96 


4.7 xlO^ 


0.992 


5 xlO"'' 


8 xlO"^ 


-2 xlO"" 


1 


-0.032 


16 


0.91 


0.03 


0.03 


0.02 



Table 2. Thermal correction to the pressure between Au 
plates at T = 300 K (column 2) and different contributions to 
it as a function of separation computed using the impedance 
of infrared optics Zp. See text for further discussion. 



a (/im) 


AP (mPa) 


AP 


^-Pte,ew 


"^-Pte.pw 


^-Ptm,ew 


^-Ptm.pw 




AP 


AP 


AP 


AP 


0.2 


-0.0097 


4.7 


-0.76 


0.10 


1.55 


0.1 


0.25 


-0.0081 


4.0 


-0.43 


0.13 


1.18 


0.12 


0.3 


-0.0070 


3.4 


-0.26 


0.14 


0.97 


0.14 


0.35 


-0.0059 


2.9 


-0.18 


0.17 


0.84 


0.17 


0.4 


-0.0051 


2.5 


-0.13 


0.20 


0.73 


0.20 


1 


-0.0026 


1.3 


-0.01 


0.39 


0.2 


0.4 



Table 3. Values of the emittivity e of Au at T = 295 K. 



Measured 


Lifshitz Z = Zn 

(e = Ed) 


Z = Zd 


Z = 


Zt 


Z — Zp 


0.02 


0.0098 0.023 


0.0098 


0.013^ 


■0.016 


2.5 xlO"* 



